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Abstract 

Full one-loop electroweak corrections to the on-shell decay H + — ► I^ + j4 are com- 
puted in the framework of models with two Higgs doublets (THDM). Such a decay 
may be dominant for over a wide range of parameter space relevant at present 
and future colliders. We show that the corrections may approach 40% and in partic- 
ular are sensitive to A5, which parametrizes the discrete symmetry breaking term. 
We suggest that a measurement of the branching ratio of H + — ► PF + ^4° may offer 
a possibility of measuring the magnitude of A5. 
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1. Introduction 



The phenomenology of charged Higgs bosons (i? ± ) has received much attention in recent 
years [[[[] since their discovery would provide conclusive evidence of physics beyond the 
Standard Model (SM) 0. Charged Higgs bosons are predicted in many theoretically 
well-motivated extensions of the SM. The simplest model which contains a if ± is the Two 
Higgs Doublet Model (THDM), which is formed by adding an extra complex SU(2)l ® 
U(l)y scalar doublet to the SM lagrangian. Motivations for such a structure include 
CP-violation in the Higgs sector and a possible solution to the cosmological domain wall 
problem ||. In particular, the Higgs sector of the Minimal Supersymmetric Standard 
Model (MSSM) takes the form of a constrained THDM. 

The phenomenology of has received substantial attention at e + e~ colliders 0,0, 
hadron colliders ||, 0, ||, colliders || and 77 colliders [[LTJ. Most phenomenologi- 

cal studies have been carried out in the context of the MSSM. The combined null-searches 
from all four CERN LEP collaborations derive the lower limit m#± > 77.4 GeV (95% c.l) 
11], a limit which applies to all models in which BR(i/ ± — > rz/ T ) + BR(/7 ± — > cs)=l, 



where BR signifies branching ratio. Current mass bounds from LEP-II for the neutral 
pseudoscalar A of the MSSM (m A > 90.5 GeV) force m H ± > 120 GeV in this model ][Lj, 
which is stronger than the direct search limit above. Limits on m#± from the Fermilab 
Tevatron searches |13| are tan j3 dependent since a significant BR(t — > H + b) is required 
in order to obtain a visible signal. The limits are competitive with those from LEP-II for 
the regions tan (3 < 1 or > 40. 

In the MSSM, m#± and the mass of the pseudoscalar are approximately degenerate 
for values greater than 200 GeV, and so the two body decay H ± — > A°W is never allowed 
for masses of interest at the CERN Large Hadron Collider (LHC). The three-body decay 
H — > A°W* — > A°ff is open for smaller m H ± although it possesses a small branching 
ratio (BR< 5% for m#± > 110 GeV) [p~4[| . Therefore for masses of interest at the LHC the 
principal decay channel is H ± — > tb, with decays to SUSY particles possibly open as m#± 
increases [15]. Recently there has been a surge of interest in studies of the MSSM with 



unconstrained CP-violating phases. In such scenarios H° — A mixing may be induced 
radiatively ][ilj although this only leads to maximum mass splittings of the order 20 GeV 
for small values of tan/?; thus the two body decay would not be open even in this case. 

Non-supersymmetric THDMs (hereafter to be called simply 'THDM') have also re- 
ceived considerable attention in the literature. In such models all the Higgs masses may 
be taken as free parameters (in contrast to the MSSM), thus allowing the possibility of 
the two body decay — > A°W for certain choices of and m^±. This decay mode 
posesses no mixing angle suppression, in contrast to — > h°W, and may compete with 



conventional decays ||17||,||18||. In fact, in a sizeable region of parameter space we show 



that it may be the dominant channel. Motivated by the results in Ref. [17] the authors 



calculated in Ref. |T9[ the Yukawa corrections to the decay H ± — > A !-!^*) and found that 



in the on-shell case the corrections may approach 50% for small values of tan (3. In this 
paper we complete that analysis and include the full bosonic corrections for the case of 
the W being on-shell. 

Conventional Higgs searches at LEP-II assume the decays H ± — > ru T and cs fly]], 
although the OPAL collaboration has recently carried out the first search for H ± — > A°W* 
topologies [^0| in the context of the THDM (Model I). A recent study [0] showed that the 
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decay — > A !^ offers good chances of detection for H at the LHC, where an analysis 
in the context of the MSSM with an extra singlet superfield was carried out (NMSSM). 
Much of this work would be relevant for the THDM that we consider, although our bosonic 
corrections would not be directly applicable to the NMSSM, since the latter possesses two 
more neutral Higgs bosons in addition to different Higgs self-couplings. 

The paper is organized as follows. In section 2 we introduce our notation and outline 
the form of the 1-loop corrections. In section 3 we explain the importance of including 
diagrams with the emission of a soft photon, H + — > W + A 0/ y, in order to keep the radiative 
corrections infra-red finite. Section 4 covers the various experimental and theoretical 
constraints that we impose. Section 5 presents our numerical results for the full corrections 
(bosonic and Yukawa), while section 6 contains our conclusions. The explicit form of the 
corrections is contained in the appendix. 

2. Lowest order result and structure of one- loop ra- 
diative corrections 



2.1 Lowest order result 



We will be using the notation and conventions of our previous work p9| , which we briefly 
review here. The momentum of the charged Higgs boson H + is denoted by pn {ph is 
incoming), pw is the momentum of the W + gauge boson and pa the momentum of the 
CP-odd A (pw and pa are outgoing). 

The relevant part of the lagrangian describing the interaction of the W with if ± and 
A comes from the covariant derivative which is given by: 

C = -£-W+(H- t A ) +h.c. (2.1) 

This interaction is model independent (SUSY or non-SUSY) and depends only on stan- 
dard parameters: electric charge (e) and Weinberg angle (s-w = sin9 w ). 

The lowest-order Feynman diagram for the two-body decay H + — > is depicted 

in the following figure: 




Figure. 1 



In the Born approximation, the decay amplitude of the charged Higgs into an on-shell 
CP-odd Higgs boson A and the gauge boson W + (Fig.l) can be written as: 

M°(H+ - W+A°) = e ;r£ where rft = %^—{p H +Pa)» (2.2) 



3 



Here e„ is the W polarization vector. We then have the following decay width: 



T °on = 1fi 2 °2 ( m ^± . m A> ™w) (2-3) 

16s w m w m 6 H± 

where A = X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) is the familiar two-body phase space 
function. Note that in the MSSM the two-body decay of the charged Higgs boson into 
W /+ v4° is kinematically not allowed. In this paper we will not present results for the case 
of W being off-shell. Ref. [19| evaluated the Yukawa corrections in the off-shell case for 



rriH± i n the range of LEP-II, finding maximum values of a few percent. 
2.2 One— Loop radiative corrections 

We shall evaluate the bosonic one-loop radiative corrections to the decay H + — > W + A°, 



and add them to the Yukawa corrections previously evaluated in Ref. [19|. This set of 



corrections is ultra-violet (UV) and infra-red (IR) divergent. The UV singularities are 



treated by dimensional regularization |22| in the on-mass-shell renormalization scheme. 
The IR divergences are treated by the introduction of a small fictitious mass S for the 
photon, which we shall explain in the next section. 

The typical Feynman diagrams for the virtual corrections of order a are listed in figure 
2.1 — >2.16. These contributions have to be supplemented by the counterterm renormal- 
izing the vertex H + A W~ (eq [21]). Note that in the THDM, the vertices W + A°G~, 
W + G°H~, W + W~A° and A°H + H~ are not present, and so the mixing G + -H + , G°-A° 
and W + -H + does not give any contribution to our process. In our case, the gauge boson 
W is on-shell and so the mixing W ± -G T is absent. The full set of Feynman diagrams 
are generated and computed using the FeynArts and FeynCalc |23], |24j packages. The 



amplitudes of the typical vertices are given in terms of the one- loop scalar functions [25 



and are written explicitly in appendix B. We also use the fortran FF-package |26| in the 
numerical analysis. 



In what follows we will use the on-shell renormalization scheme developed in Ref. |L9 
(and refs therein). The vertex counterterm is given by: 

8C = -^W+(H~ 7 A°)(~5Z WW + hz A o AO + Uz H±H± + 6Z e - ^) (2.4) 

ZS\y Z Z Z S\y 

where 5Zy/w, ^Z^oa and 5Zh±h ± are the wave function renormalization constants for 
the W ± gauge boson, A and H ± Higgs boson defined as follows: 

8Z tt = -^lP-\ k2=mf z = W,A° H* (2.5) 
Sm 2 = ReY> u (m 2 ) i = W, Z (2.6) 

where Ejj(fc 2 ) is the bare self-energy of the H ^ , A or W. The electric charge counterterm 
and are defined as: 

SW 



(2.7) 



Ss w _ 1 c 2 w bm 2 w _ 5m% 

r> 2 V 2 1 I 

S W ZS W 77% 777f 

A7 K 7 ls w I dYPpjk 2 ) , s w gg(0) r _, 

5Z e = --8Z,, + -— 5Z Zl = 2^T^U 2 =o + — {2.1 
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The index T in Sj 7 and E 7 -^ denotes that we take the transverse part. The scalar one- loop 
self-energies entering in the above equations » pT8|) are given in appendix C\ and C 2 , 
while gauge boson self energies can be found in [ p^jj . 

The one- loop amplitude A4 1 (vertex plus counterterms) can be written as: 

M \H+ - W + A°) = ^{Thp'h + T WI f w )e; (2.9) 

where Th and Tw can be cast as follow: 

IV = T v ^ rtex + 5T v ^ rtex (2.10) 
T H = F™ rtex + 5T v * rtex (2.11) 

Here T^ r ^ x represents the vertex corrections and ST^ r ^ x is the counterterm contribution 
needed to remove the UV divergences contained in T^ r jf x . 
The expressions of the counterterms are: 

6v verte X = _ S JW + h 5Zu+H+ + ^ + ^)) 

^vertex = _ 2 ^ HeX (2.12) 

In the on-shell case the interference term 2ReA4°*Ai 1 , found from squaring the one- loop 
corrected amplitude \A4° + A4 1 \ 2 , is equal to r#|.M | 2 fL9"f . Hence the one-loop corrected 
width Tl n can be written as 

T 1 on = (l + T H )T° on (2.13) 

with r# (defined by eq. |2.11|) being interpreted as the fractional contribution to the 
tree-level width, T® n . Note that IV (eq. 2.10) does not contribute to Tj n . 



3. Real photon emission: H + —> W + A°^f 

The vertex correction supplemented by the counterterms is UV finite but there still re- 
mains infra-red divergences. These arise from the diagrams 2.10 and 2.11 with V = 7 
and also from the wave function renormalisation constant 5Zh±h ± (Diagram 2.33) and 
5Zy/w In order to obtain a finite result, one has to add the correction from the emission 
of a real photon in the final state as drawn in figures 2.17 — > 2.19. 

In terms of the momenta of the particles in the final state, the square amplitude of 
the process H + — > A Q W +r y is given by: 

|M(tf+^A W/+ 7 )| 2 = ~\{m\ ± ,m\^)[^ E± 1 



Syy Ai ' W/L jyi^ {2pHki) 2 



' {2p w k 1 ) 2 ' v m 2 w l \2p w k 1 ){2p H k 1 ) [ (3-1) 

where fc 7 denotes the momentum of the photon. Note that as a consequence of gauge 
invariance, the amplitude of the sum of the three diagrams (Figs. 2.17 — > 2.19), should 
vanish when multiplied by the four-momentum of the photon, which provides a good 
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check of the calculation. The integrals over three body phase space can be found in [25], 
and one obtains the following expression for the width: 



o e 2 M 2 H± m 2 H± rn 2 H ±-m\ 

J- Br — — 1 on 77J T~ 1 HH + IWW + I- 1 H o )J-HW\ (A^J 

tc 2 \2 m w mfo 
Where Ihh , Iww an d Irw are given as follows: 

1 ril, / ^ \ \1 / 2 2m /A 



W = -^{2\og(- - )log(/3 2 ) + 21og 2 (&) - log 2 (/3 ) - log 2 (A) 

Am H± om H ±mAm w 

+2Sp(l - Pi) - Sp(l - /3 2 ) - Sp(l - A 2 )} (3.3) 

where A = \(m H ±, m\, rriyy) is the two body phase space, 5 is a small fictitious photon 
mass, Sp is the dilogarithm function and /3, are defined as: 

_ m 2 H± - m 2 ^ - + m^± - m 2 w + m 2 A - 

Po ~ 2^r A ' Pl ~ ' 



j3 2 = H± + m W~ m A~^ /g 4 x 

2m H ±m w 

We stress here that the IR divergence contained in Ihh (Iww) is cancelled by the wave 
function renormalisation constant of the charged Higgs H ± (W^, while the IR divergence 
contained in Ihw is cancelled by the vertex diagrams 2.10 and 2.11 (with V = 7). One can 
confirm easily that adding the virtual corrections with the Bremsstrahlung diagrams yields 
an IR finite result. This feature has been checked both algebraically and numerically. 

4. THDM scalar potential: Theoretical and Experi- 
mental constraints 

In this section we define the THDM scalar potential that we will be using. In appendix A 
we list the trilinear and quartic scalar self-couplings which are relevant to our study. Other 
relevant couplings involving Higgs boson interactions with gauge bosons and fermions 
can be found in Ref. J5J. For a full list of scalar trilinear and quartic couplings see 
Ref. [F2UJ. It has been shown |TJ that the most general THDM scalar potential which is 
both SU(2)l <8> U(l)y and CP invariant is given by: 

Vfr, $2) = \i{\*i\ 2 - vl) 2 + A 2 (|$ 2 | 2 - v\f + A 3 ((|$!| 2 - vl) + (|$ 2 | 2 - v 2 2 )) 2 + 

A 4 (|$i| 2 |$ 2 | 2 - \^2?) + A 5 (i?e($+$ 2 ) - v lV2 f + A 6 [im($+$ 2 )] 2 (4.1) 

where $1 and $ 2 have weak hypercharge Y=l, v± and t> 2 are respectively the vacuum 
expectation values of $1 and $ 2 and the A, are real-valued parameters. Note that this 
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potential violates the discrete symmetry $j — > — $j softly by the dimension two term 
\z,Re(&i&2) and has the same general structure as the scalar potential of the MSSM. 
One can prove easily that for A5 = the exact symmetry $j — > — $j is recovered. We 
note that Ref. 0] lists the complete Higgs trilinear and quartic interactions for two 
6 parameter potentials, referred to as 'Potential A' and 'Potential B'. Potential A is 
equivalent to our potential if A5 — > 0, and in this limit the Feynman rules in the appendix 
A are in agreement with those in Ref. pi . 



After electroweak symmetry breaking, the W and Z gauge bosons acquire masses given 
by m w = \9 2v2 an d fn\ = \i.9 2 + fi ,/2 ) t ' 2 ? where g and g' are the SU{2) L and U(l)y gauge 
couplings and v 2 = v\ + v\. The combination v\ + v\ is thus fixed by the electroweak scale 
through v\ + v\ = (2^/2Gp)~ 1 , and we are left with 7 free parameters in eq.(4.1), namely 



the (Aj)^!,...^ and tan/3 = i^/Vl- Meanwhile, three of the eight degrees of freedom of the 
two Higgs doublets correspond to the 3 Goldstone bosons (G ± , G°) and the remaining five 
become physical Higgs bosons: if , h° (CP-even), A (CP-odd) and H ± . Their masses 
are obtained as usual by the shift $. ; — > $j + v,- t and read ffl]: 



m 2 A = X 6 v 2 , m 2 H± = X 4 v 2 , m 2 H h = ±[A + C ± sj '{A - C) 2 + 

where 

A = 4u 2 (Ai + A 3 ) + v 2 X 5 , B = ^ 2 (4A 3 + A 5 ) and C = Av 2 2 {\ 2 + A 3 ) + ^ 2 A 5 (4.2) 

The angle f3 diagonalizes both the CP-odd and charged scalar mass matrices, leading to 
the physical states H and A . The CP-even mass matrix is diagonalized by the angle 
a, leading to the physical states H°, h°, with a given by: 

2B A-C , . 

sin 2a = , = , cos 2a = = (4.3) 

yJ(A - C) 2 + AB 2 ,J(A-C) 2 + 4B 2 

It is then straightforward algebra to invert the previous equations to obtain the Aj in 
terms of physical scalar masses, tan/3, a and A5: 



\ 9 2 2 \ 9 2 2 \ 9 2 s aC a , 2 2\ A5 

A4 = 2^ m ^ ' K = 2^ mA ' A3 = 8<:^c; (m --^ ) - t (44) 

Ai = ld-rKm 2 H + slml - ^-(m 2 H - m 2 h )} - ^(-1 + tan 2 /?) (4.5) 



n 2 \ 

y r_2 2 , _2„_2 _ _ /3/„_2 „_2m a 



A2= qJ^2 i s i m H + cl m h ~ s a c a tan f3(m 2 H - m{)} - ^(-1 + — ^— ) (4.6) 
oSpTn-^y q tan p 

We are free to take as 7 independent parameters (Aj)j = i r .. i6 and tan/? or equivalently the 
four scalar masses, tan/3, a and one of the A,. In what follows we will take A5 as a free 
parameter. In our analysis we also take into account the following constraints when the 
independent parameters are varied. 

• The contributions to the 5p parameter from the Higgs scalars BIJ] should not exceed 



the current limits from precision measurements [32]: —0.0017 < 5p < 0.0027. 



From the requirement of perturbativity for the top and bottom Yukawa couplings . 
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tan/3 is constrained to lie in the range 0.3 < tan (3 < 130. Upper and lower bounds have 
also been obtained from the experimental limits on the processes e + e~ — > Z* — > h 0/ ~f 
and/or e + e~ — > A°^. For very light h or A (« 10 GeV) Ref. |!4| derived 0.15 < tan/3 < 
75, with the limits weakening for heavier m/^m^). For our study we will restrict the 
discussion to values tan/3 > 0.5. 

• We require that tree-level unitarity is not violated in a variety of Higgs scattering 
processes |BH . 



5. Numerical results and discussion 

In this section we present our numerical results for Th, which is the fractional correc- 
tion to the tree-level width (eq.2.13). We take the following experimental input for the 
physical parameters [P^j . The fine structure constant: a = f- = 1/137.03598, the gauge 
boson masses: mz = 91.187 GeV, mw = 80.41 GeV, the lepton masses: m e = 0.511MeV, 
m M = 0.1057GeV, m T = 1.784GeV. For the light quark masses we use the effective val- 
ues which are chosen in such a way that the experimentally extracted hadronic part of 
the vacuum polarizations is reproduced [RTJ: m d = 47 MeV , m u = 47 MeV , m s = 
150 MeV , m c = 1.55 GeV , = 4.5 GeV. For the top quark mass we take m t = 175 

2 

GeV. In the on-shell scheme we consider, sin 2 9w is given by sin 2 6\y = 1 — r ^-, and this 
expression is valid beyond tree-level. 



Let us make some comment about the Yukawa corrections discussed in Ref. |TjJ. In 
the case of an on-shell W , we showed that for small tan /3 in both Model I and II one can 
find large corrections of up to around ±10% (away from threshold effects at ttla ~ 2m t ). 
In Model I for large tan/3 all fermion corrections decouple and reach a constant value 
of 3.3% for tan/? > 4. In Model II Th is enhanced for large tan/3 (> 20) since in this 
scenario the internal b quarks in the loop couple more strongly to A . Typically for 
rriH± = 440 one can reach a correction of about —7% — > —20% for large tan/3 > 90 and 
light irtA (60 < iriA — 100 GeV). For > 100 GeV the tan/3 dependence of the Yukawa 
correction is rather weak and lies in the range of —5% — > 5% for 100 < rriA < 330 GeV. 

We stress at this stage that in the low tan/3 (tan/3 < 2) regime the corrections in 
Model I and Model II are practically identical. Perturbative constraints on the A« and 
unitarity constraints on the quartic scalar couplings constrain the magnitude of tan /3. As 
shown in |35|, tan/3 > 20 violates the unitarity bounds if A 5 = 0, although for A 5 ^ 



values of tan/3 > 40 are comfortably allowed. However, perturbative constraints on the 



Aj (in particular Ai) disfavour tan/3 > 30 | 3~5| , and so in our analysis we will only consider 



small to moderate values of tan/3. Therefore our results are applicable to both Model I 
and II. Note that in order to satisfy the experimental constraint on 5p we have assumed 
(for the graphs we have plotted) that a — /3 — | and the charged Higgs boson mass m H ± 
is quasi-degenerate with mu- 

In Fig. 3 we plot Y H as a function of m H ± for m H = m H ± — 10 GeV, rrih = 120 GeV, 
rriA = 150 GeV and a = /3 — | for several values of A 5 . With the above set of parameters 
and for tan/? = 0.5(1.5) and A 5 = in Fig. 3. a (Fig.3.b), the unitarity constraints in 
the spirit of Ref. |35[ require m#± < 370 (480) GeV, with this bound weakening for 
increasing A5. This can be seen both in fig.3.a and fig.3.b, where we cut the curves at 
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the value of the charged Higgs mass which violates the unitarity constraint. Both for 
Fig.3.a and Fig.3.b the Yukawa correction is positive and lies in the range 3.6% — > 11.9% 
and 3.3% — > 4.3% for Fig.3.a and Fig.3.b respectively while the bosonic correction is 
negative. In Fig.3.a (tan/3 = 0.5), the bosonic correction is in the range —0.6% — > —0.9% 
for A5 = and m H ± G [231, 370], and so in this case the Yukawa correction is dominant; 
for larger A 5 = 8.5 the bosonic correction becomes strongly negative ( —11% — > —50% for 
m H ± G [230,620]) and dominates the Yukawa corrections. 

In Fig.3.b we take tan/3 = 1.5, and the bosonic correction is in the range —0.6% — > 

— 1.6% for A 5 = and m H ± G [230, 475]; for A 5 = 1 and m H ± G [230, 520] it is in the range 

— 1.7% — > —4.7% . In the above two cases the bosonic and Yukawa corrections interfere 
destructively leading to a small total correction of about pa 3%. In the case where A 5 > 3 
the bosonic correction becomes strongly negative and dominates the Yukawa correction. 

Fig. 4. a and Fig.4.b show the total contribution to Th as function of itla for mn = 
500, rrih = 360 and m H ± = 530 GeV. Fig. 4. a corresponds to tan/3 = 0.8 and Fig.4.b 
corresponds to tan/3 = 1.6 (with a — f3 — f). In both figures the Yukawa correction is 
positive in the region niA < 277 GeV and itia > 352 GeV, while for rriA ~ 350 GeV the 
channel A — > tt opens, leading to a very large negative correction. Note that the bosonic 
correction is negative for every value of uia- 

We can conclude that for the intermediate mass range of itla (away from threshold 
effects niA ~ 2m t ) and for A 5 < 4 there is a cancellation between the Yukawa correction 
and the bosonic correction. For large A 5 the contribution to T H is dominated by the 
bosonic correction and is consequently negative. For itla ~ 350 GeV (~ 2mt) there is 
constructive interference between the Yukawa and bosonic corrections. 

In Fig. 5. a we plot the bosonic contribution to Th (denoted T b ^ s ) as a function of A5 for 
tan/3 = 0.5,5,8. The fermionic correction (independent of A 5 ) takes the values -3.21%, 
3.22%, 3.24%. One can see from the curves that T b ^ s increases with A 5 and may approach 
—40%. This is expected from the form of the trilinear scalar couplings which increase 
linearly with A5. 

Fig.5.b shows the dependence of T^ 5 on sin a for — n/2 < a < ir/2. Here C±, Ci and 
C 3 correspond to three distinct parameter configurations which satisfy both the unitarity 
and p parameter constraints (see figure caption). The fermionic corrections (independent 
of a) for Ci, C 2 and C 3 take the values 3.2% , 3.5% and 3.4% respectively. The bosonic 
corrections become important for sin a pa ±1. 

It is apparent from the figures that r^ s has a complicated dependence on a, /3, A5 and 
the physical Higgs masses. This is clear from the explicit form of the trilinear couplings 
in appendix A.l, which mediate the numerous triangular loop corrections. Therefore 
enhancement in the bosonic sector may occur in a variety of scenarios. Of particular 
interest is the sensitivity to A5, a measurement of which (along with the Higgs masses 
and mixing angles) would allow reconstruction of the Higgs potential. We suggest the 
measurement of BR(i7 ± — > A !^) as a way of obtaining information on A 5 . The decay 
mode {H^ — > A°H^) may in fact be dominant and in Fig.6 we show the ratio 

R _ Eg - A ° W ^ (5 1) 

r(#± -> a°w) + r(#± -> tb) 1 ; 

as a function of tan/3 for various values of itla, fixing m^± = 500 GeV. We plot the tree- 
level width for (H^ 1 — > A°W) given in Eq. 2.3, and the tree-level width for H ± — > tb given 
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in Ref. [T%\, assuming Model II type couplings. Since other channels such as H ± — > h°W 
and H ± — > H°W may be open the above ratio should be interpreted as the upper bound 
on BR(i/ ± — > A°W). Note that these additional channels are suppressed by the factors 
cos 2 (/3 — a) and sin 2 (/3 — a) respectively, while r(if ± — ► possesses no mixing 

angle suppression. One can see that the decay — > A°W is maximized for moderate 
values of tan/3 (i.e. when r(if ± — > tb) is minimized). The curves with lighter 
are less phase space suppressed and so the value of R may be larger. For fixed m A 
and tan/? the sensitivity to m H ± is rather mild. Larger m#± slightly increases R since 
r( J ff ± -> A°W) ~ m 3 H± while T( J ff ± -> t6) ~ m H ±. 

Given the possible large BR, an accurate measurement of BR(ff ± — > A ! 7 ^) may allow 
one to obtain information on A5. As explained above, the radiative corrections show 
sensitivity to several of the input parameters. If experimental information on the Higgs 
masses and mixing angles were available then it might be possible to measure A5. At a 
e + e~ linear collider [^HJ one could measure the Higgs masses from a variety of production 
mechanisms (e + e~ — > Zh° , ZH°, h°A° , H°A° , H + H~). Information on the mixing angles 
a, (3 could be obtained |5I5| from an analysis of the production cross-sections and branching 
ratios. Other processes which are sensitive to A5 are e + e~ — > H + H~ || and H ± W T [[40f 



while theoretical bounds on the Higgs masses in the case of A5 7^ are explored in 
Ref. pLH. 



6. Conclusions 

We have computed the radiative corrections to the on-shell decay H + — ► A ^" 1 " in the 
general Two Higgs Doublet Model, taking into account the experimental constraint on 
the p parameter and also unitarity constraints on the scalar sector parameters. We have 
included the Yukawa corrections, the full electroweak corrections (bosonic), and also the 
real photon emission in the final state (Bremsstrahlung). The computation was done with 
dimensional regularization in the on-shell scheme. We find that the total radiative correc- 
tions may approach 40% in regions of parameter space for both small and moderate tan/3. 
The bosonic correction is sensitive to the soft discrete symmetry breaking parameter A5, 
and may interfere both constructively and destructively with the Yukawa correction. For 
larger A5 the bosonic contribution becomes strongly negative and in general dominates 
the Yukawa correction. For tua ~ 2m t and low tan (3 the Yukawa correction is maximized 
and interferes constructively with the bosonic correction, resulting in large negative cor- 
rections to the tree- level width. Finally, we showed that the decay H ± — > A ^ 1 * 1 may 
supercede H — > tb as the dominant decay channel, and thus a precise measurement of 
its branching ratio may allow information to be obtained on A5. 



Acknowledgements 

A.G.A was supported by the Japan Society for Promotion of Science (JSPS). We thank 
C. Dove for reading the manuscript. 



10 



Appendix A: THDM trilinear and quartic scalar cou- 
plings 

In this appendix we list the Feynman rules in the general THDM for the trilinear and 
quartic scalar couplings relevant for our study. All formulae are written in terms of the 
physical masses, a, f3 and the soft breaking term A 5 . Note that in the trilinear couplings 
(quartic couplings) we have factorised out ie (ie 2 ). In the following g c = 1/ {2swmw s 2p) 
, v 2 = , Cp a = cos((3±a) and s^ a = sm(/3±a). 

A.l Trilinear scalar coupling 



9hoh+h- = -2g c (m 2 H o(c 3 pS a + s^c a ) + m 2 H± s 2fi Cp a - s^ a X 5 v 2 ) (A.l) 

9h°h+g- = 9cS2pSp a (m 2 H0 - m 2 H± ) (A.2) 

9h°H+H- = -2g c (m 2 h0 (c a cp - s a s 3 p ) + m 2 H ±s 2 pSp a - c^ a \ 5 v 2 ) (A.3) 

9h°H+G- = -9cS2f3Cp a (m 2 h o -m 2 H ±) (A.4) 

9h°a°a° = -2g c (m 2 H0 (s a c 3 p + c a sj|) + m 2 A s 2(3 c^ a - s^ a \ 5 v 2 ) (A.5) 

9h°a°g° = 9cS2pSp a (m 2 H0 - m 2 A ) (A.6) 

9h°A«A° = -2g c (m 2 h0 (c a cl - s a sj) + m A s 2p-Sp a - c t^v 2 ) (A.7) 

gh°A°G° = 9cS2/3Cp a (m 2 A - m 2 h0 ) (A. 8) 

9a°h+g- = ig c S2p(m 2 H ± - m 2 A ) (A.9) 



A.2 Quartic scalar coupling 



9h°h°h°h° = -l2g 2 c [m 2 H0 {c p sl + spcl) 2 + m 2 h0 s 2 a c 2 a s p 2 - A 5 v 2 (c 2 - cj) 2 ] (A.10) 

gh°h°h°h° = -l2g 2 c [m 2 H0 s 2 a c 2 a c^ a 2 + m^cgc* - s p s 3 a ) 2 - \ 5 v 2 (c 2 a - s 2 ^ 2 } (A.ll) 

9a a°a°a = -l2g 2 c [m 2 HO (c a sl + s a cj) 2 + m^ (c a c| - s a Sp) 2 - \ 5 v 2 c 2 2 p} (A. 12) 

g G °G°G°Go = -3g 2 c s 2 2 p[m 2 H oCp a 2 + m 2 h0 sp 2 ] (A.13) 

9H+H-H+H- = -85-c[ m io( C a4 + S«4) 2 + m M»( C « C /3 _ S "4) 2 ~ A 5^ 2 C^] (A. 14) 

9G+G-G+G- = -2g 2 c s 2 w [m 2 HoC - p 2 + m^ s^ 2 ] (A.15) 
gw>H"h°h° = -4gc[m 2 H0 s a c a (spCf3 + 3s a c a Sp a 2 ) - m 2 h0 s a c a (spcp - Ss a c a c^ a 2 ) 

-X 5 v 2 (3s 2 a c 2 a -s 2 c 2 )} (A.16) 
gH H°A°A = -2g 2 c [2m 2 H0 (s a cl + c a s z p ){cps z a + spc 3 a ) + 2m 2 h0 s a c a Sp a (c a c 3 p - s a sj) 

+m 2 4 (cj - sl)(c 2 a -sl + 2s a spcp a ) 

+2A 5 t; 2 (c 2/3 (c 2 4 - s 2 a 4) - S y a (l + s 2a s 2 p))] (A.17) 
gh°h°AOA° = -2g 2 c [2m 2 H0 {s a c ? p + c a s 3 p )s a c a Cp a + 2m 2 hQ (c p c ? a - sps 3 a )(c a c 3 p - s a sj) 
+m\{<?p - s 2 p ){s 2 a + s\ - 2s a s p cp a ) 

+2A 5 t; 2 (c 2/3 ( S 2 4 - c 2 a 4) - S y a (l - s 2a s 2 p))] (A.18) 
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Appendix B: One-loop vertex: H + — > WA° 

In this appendix we will list the analytic expression for each generic diagram of Fig.2. The 
sum over all the particle contents yields the corresponding contribution to the one-loop 
amplitude for H + — > WA°. For scalar and tensor integrals we use the same convention 
as P, the analytical expression of all scalar functions can be found in P5L E?J. 

B. 1 Fermionic Loops 

The diagram with the uud triangle, Fig. 2.1, yields the following contribution to the one- 
loop amplitude: 

M2.1 = - 2 e jf C sw Y uu {Y u L d (2B + B x ) + m u C (m u Y u L d + m d Y u R d ) 

+Y u L d [m 2 w C 2 - 2C 00 - 2m 2 A C n + C 12 (m 2 H± - m 2 w - m 2 A )]} (B.l) 

Nc = 3(1) for quarks (leptons). Therein, all the B , Bi, Ci and Cy have the same 
arguments: 

[B^B^ml^mlm 2 ) , [C h dj] (m 2 A , m 2 H ± , p 2 w , m 2 u , m 2 u , m 2 d ) 

The summation has to be performed over all fermion families; in practice only the third 
quark generation is relevant. 

The corresponding expression for the the diagram with the ddu triangle in Fig. 2. 2 is 
obtained from the previous one by making the following replacements: 

Y ud < ► Y ud > Y uu < ► Y dd , m u < ► m d 

B. 2 Bosonic Loops 
Diagram 2.3 

The couplings are summarized in the following table: 
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( Q Q Q \ 

yiJii Jji Ok ) 


9H+SiS k 


9ASiSi 


9w+s,s k 


[n, A, a ) 


9hH+H+ 


9hAA 


— zfi — 


(G+,H+,h) 


9hH+G+ 


9AG+H+ 


2s w 


(G + ,H+,H) 


9HH+G+ 


—igAH+G+ 


S l3a 

2s w 


{G+,H+,A) 


—WAH+G+ 


—i9AH+G+ 


I 

2sw 


(HA lT+\ 


9HH+H+ 


9HAA 


1 

2sw 


{h,G°,G + ) 


9hH+G+ 


9hAG° 


1 

ZS\Y 


(if, G°, G+) 


9HH+G+ 


9hag° 


1 

2s\.\r 


(A,h,G+) 


—WAH+G+ 


9hAA 


pot 


(A, H, G+) 


WAH+G+ 


9haa 


pa 

2sw 


(H+,G+,h) 


9hH+H+ 


—i9AH+G+ 


s7 

pa 

2sw 


(H+,G+,H) 


9HH+H+ 


—i9AH+G+ 


pot 
2s\y 



Diagram 2.4 



M 2A [h] = -^^{-±B {m 2 H± ,m 2 h ,m 2 w ) + [2{-m 2 A + m 2 H± -2m 2 z )- 4p* w ]C [h] 

+[3(-m 2 A + m 2 H± )- Ap^C^h] + [A(-m 2 A + m 2 H± ) - 4p 2 w }C 2 [h] + 4C 00 [h] 
+ [m 2 H± + 3m 2 A - p 2 w ]C n [h] - 2[m 2 H± - m 2 A ]C 12 [h}} (B.3) 
M 2A [H\ = M 2A [h][cp a 2 ^Sp a \m h ^m H ,h^H] (B.4) 



Where the arguments of Cj and are given as follows: 

[Ci,Cij}[S] = \C i ,C i3 \{m 2 A ,m 2 H± ,p 2 w ,m 2 z ,m 2 s ,m 2 w ) 



Diagram. 2. 5 



M 2 . 5 [h] = -— -^{B (m 2 H ±,m 2 w ,m 2 h ) + Bi(m^±, m^,m^) + [p^ - 2m 2 H ±]C 1 [h] 

L07TS W 

+2C 00 [h] + [m 2 H± - (p 2 w - m 2 A )]C n [h] - [Vw -™ 2 A + m 2 H± ]C 12 [h]} (B.5) 
M 2 . 5 [H] = M 2 . 5 [h][cp a ^ Sp a ,m h -+m H ,h^ H] (B.6) 

M 2 . 5 [A] = M 2 4h][cp a 2 ^-l,m h ^m A ,h^A] (B.7) 
Where the arguments of Cj and Cij are as follows: 

[Q, C^] [S] = [d, Cij] (m A , m 2 H± , p 2 w , m 2 H ± , m 2 w , m 2 s ) 
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Diagram. 2. 6 

87TC^ 

M 2 . 6 [tf] = M 2 . 6 [h][cp a ^-Sp a ,h^H] (B.9) 
Where C and C\ have the same arguments as follows: 

[Co , Ci][S] = [C , Ci](m^,m^±,p^,m|,m|,m^) 

Diagram. 2. 7 



6Ct(c^ )c _ 

+[& - ™h - "Ih±]Ci[/i] + 2C 00 [/i] + [-Pw + m\ + m 2 H± ]C n [h] 
-\p 2 w -m 2 A + m 2 H ±]C 12 [h}} (B.10) 
M 2 . 7 [H] = M 2 . 7 [h][cp a ^ sp a ,h-^ H,m h ^m H ] (B.ll) 

The arguments of Cj and Cy are given by : 

[CudjWS] = \C h C i3 \{m\m 2 H± ^ 2 w ,m 2 s ,m 2 z ,m 2 H± ) 

Diagram. 2. 8 

M 2 . 8 [/i] = eaS f a ~ 2 mw g h o AOA o(2C [h] + CM) (B.12) 
M 2 . 8 [iJ] = M 2 . 8 [/i] [s^ Q -> c^ Q , H, g h o A o A o -> s^o^o] (B.13) 
Where C and Ci have the same arguments: 

[C ,Ci][S] = [Co,Ci](m^,m^ ± ,p^,m|,m^,m^) 



Diagram. 2. 9 

M 2 . 9 [/i] = f^^l ffh o H+H -(2Co[/i]+C 1 [/ l ]) (B.14) 
M 2 . 9 [iJ] = M 2 . 9 [/i][s^ Q -> c^ Q , /i -> H,g h o H+H - -> g H »H+H-\ (B.15) 
Where Co and Ci have the same arguments: 

[C ,Ci][S] = [C ,Ci](m^,m^ ± ,p^,m^,m^ ± ,m|) 
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Diagram. 2. 10 



M 2 . W [Z] 



ea{(? w - s 2 w )m w 
%t\c 2 w 



^otf +G -(2C [Z] + Ci[Z]) 



M 2 . 10 [ 7 ] = £ ^^ + G-(2C [5] + C 1 [5]) 
Again the C and Ci have the same arguments: 



(B.16) 
(B.17) 



[Co , Ci][V] = [C , C 1 ](m%m%±,p^,ml [ ,,m 2 H ±,mv) 

Here 5 is a small photon mass introduced to regularise the infrared divergence contained 
in the C function. 



Diagram. 2. 11 



M 2 . 11 [Z] 



ea(c 



w 



s 2 w) 



167TS 3 



(4B (m 2 H± ,m 2 z , m 2 H± ) + [2p 2 w + m 2 A - m 2 H± + 2m 2 w }2C [Z} 



w 



+ [Ap 2 w + 3(m 2 A - m^ ± )]d[Z] + \p 2 w + m\- m 2 H± )]AC 2 [Z] - AC 00 [Z] 
+\Pw - ™ 2 H ± ~ 3rn 2 A ]Cn[Z] + [m 2 H± - m 2 A ]2C 12 [Z]) (B.18) 



M 2 . n [ 7 ] = M 2 . n [Z][ {c2w 2 -» 2, m z -> J] 



(B.19) 



All Cj and Cj,- have the same arguments 



[Ci , Cij-][V] = [C , Cij](m^,7n^±,^.,mwr,mH±,Jr4) 



Diagram. 2. 12 

M 2 .i 2 [/i,Z] 
M 2 . 12 [#+,W + ] 

M 2 . 12 [#, Z] 
Diagram. 2. 13 



eac 



_ 2 



167is w cw 



(2B (m 2 A , m 2 z , m 2 ) + ^(m^, m|, m*)) (B.20) 



167TS 



-(2B (m 2 A , m 2 w , m 2 H± ) + B x {m\, m 2 w , m 2 H± )) (B.21) 



il- 



ea s 



/3a 



167TSvi/C^ 



;2S (m^m|,m^) +S 1 (m^,m|,m|)) (B.22) 



M 2 . 13 [Z] 

M 2 . 13 [ 7 ] : 



eac 



(B (m 2 H ±,m 2 H ±,m 2 z ) - 5i(m|±, m|±, m|)) (B.23) 



87TSVK 



[B (m 2 H± ,m 2 H± ,5) - £i(m#±, m 2 H ±, 6)) 



(B.24) 
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Diagram. 2. 14 



M 2 .u = -77 — r (2B (mH±,m^,m\) + Bx(mH±,m^,m 2 A )) (B.25) 
Diagram 2.15 and 2.16 

For this kind of topology, it is clear that the amplitude is proportional to the W gauge 
boson momentum (Lorentz invariance) and consequently for the W on-shell the amplitude 
vanishes. 

M 2 . 15 = , M 2 . 16 = (B.26) 

Appendix C: Charged and CP-odd Higgs bosons self- 
energies 

This appendix is devoted to the self-energies of the charged and CP-odd Higgs bosons 
which are needed for the on-shell renormalisation scheme. The gauge bosons self energies 
7-7, 7-Z, Z-Z and W-W can be found in p7 . 



C.l CP-odd Higgs boson self-energy 

The CP-odd Higgs self-energy T, AA can be cast into three parts: i) fermionic part 2.20, 
ii) pure scalar part 2.21; 2.22, 2.23 and 2.26 iii) mixing of gauge boson and scalar 2.24, 
2.25 , 2.27 and 2.28 in such a way that: 



E^V) = (g 2 ) + 5#V) + £vs(^) (CI) 



with 



Sf (q 2 ) = ^YfYAoim^ + q'B^^m 2 )) (C.2) 

Z AA {q 2 ) = -^{g 2 hAA B {q 2 ,m 2 A ,m 2 h ) + g 2 HAA B Q {q 2 ,m 2 A ,m 2 H ) 
,2 



+9hAG B o(l ,rn h ,m z ) + g HAG B (q ,m H ,m z ) 



a 



+ 2 9ah+g- b o(c[ , m H ±,m w )) + —(-gAAAAA (m A ) - g h hAAA (m 
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-gHHAAA (m 2 H ) - 2g H+H - AA A (m 2 H± ) - 2g G+G - AA A (m 2 v ) (C.3) 

2 4 2 1 

~gAAGGA (m 2 z ) + 2 2 A (m 2 z ) + — A Q (m 2 w ) - —m 2 w - 2 2 m%) 

s w c w s w s w s w c w 

Z A j(q 2 ) = — ° 2 2 (cp a 2 (A (m 2 z ) + (m 2 h + q 2 )B (q 2 , m 2 hJ m|) 

-2q 2 B x {q 2 , m 2 , m|)) + s^ a \A (m 2 z ) + (m 2 H + q 2 )B (q 2 , m 2 H , m|) 
-2q 2 B 1 (q 2 ,m 2 H ,m 2 z )) + c^(A (m^) + (m 2 H ± + q 2 )B (q 2 \m 2 H± ,m 2 w ) 
-2q 2 BM^\^w))) (C4) 
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C.2 Charged Higgs boson self-energy 

The charged Higgs boson self-energy H H+H can be cast into three parts: i) fermionic 
part 2.29, ii) pure scalar part 2.30; 2.31, 2.32 and 2.35 hi) mixing of gauge boson and 
scalar 2.33, 2.34 and 2.36 in such a way that: 

Z H+H - (q 2 ) = Ef H± (q 2 ) + Ef *~ (q 2 ) + E^" (g 2 ) (C.5) 



^ H \q 2 ) = -^((Y^/ + Y R2 )(A (m 2 ) + q 2 B 1 (q 2 ,m 

+ (m 2 f ,(Y f L f , 2 + Y f R /) + 2m f m f Y f L fl Y f R fl )B (q 2 ,m 2 fl ,m 2 f )) (C.6) 

^S ±H± ((1 2 ) = ^9AH+G- B ^ 2 i m \i m w) + 9lH+H- B o(<l 2 , m L m H±) 

+9 2 hH+G- B o{q 2 , m 2 h , m 2 w ) + g 2 HH+H -B (q 2 , m 2 H , m 2 H± ) 
+9 2 HH+G- B o(.<l 2 , r n 2 H ,m 2 v )) - —(g H +H-AAA (m 2 A ) + g H +H-hhA {m 2 h ) 

2 

+9H+H-HHA ('m 2 H ) + 2g H+H - H+H -A {m 2 H ±) + — {m 2 w - 2A^m 2 w )) 

s w 

+9H + H-GGA (m 2 z ) + {s2w ~f ] \ m 2 z - 2A (m|))) (C.7) 



Cu/S 



" (<f) = - T ^(As 2 w ((m 2 H± +q 2 )B (q 2 ,0,m 2 H± )-2q 2 B 1 (q 2 ,m 2 H± ,0)) 



-'VS 



I6ns w 



+ (4_40! (Ao(m | ) + {m 2 H± + f)B tf, 

c w 

-2q 2 B 1 {q 2 , m 2 H± , m|)) + {A {m 2 A ) + {m 2 w + 4q 2 )B {q 2 , m 2 A , m 2 w ) 
+4q 2 B 1 {q 2 , m 2 w , m 2 A )) + c^J {A {m 2 h ) + {m 2 w + 4q 2 )B {q 2 , m\, m 2 w ) 
+Aq 2 B 1 (q 2 ,m 2 w ,m 2 h )) + s^ a (A Q (m 2 H ) + (m 2 w + 4q 2 )B Q (q 2 , 
+Aq 2 B 1 (q 2 ,m 2 w ,m 2 H ))) (C.8) 
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Figure Captions 

Fig. 1 Lowest-order Feynman diagram for the decay H + — > v4 iy + . 

Fig. 2 Feynman diagrams for the one- loop corrections to the decay H + — ► A°iy + : i) vertex 
(2.1 — > 2.16), ii) Bremsstrahlung diagrams for H + A°W + 'j: Fig. 2.17 — > 2.19 iii) 
CP-odd Higgs boson self-energy (2.20 — > 2.28) and Charged Higgs boson self-energy 
(2.29 -> 2.36). 

Fig. 3 Total contribution to as function of m H ±. We chose: m H = m H ± — 10, m h = 120, 
m A = 150 (GeV), a = (3 - § . 

Fig.3.a: tan/? = 0.5, for four values of As=0.0, 2.0, 6.0 and 8.5. 
Fig.3.b: tan/? = 1.5, for four values of A 5 =0.0, 1.0, 3.0 and 5.0 

Fig. 4 Total contribution to Th as function of tua- We chose: m# = 500, rrih = 360, 
m,H± = 530 (GeV) and a — /3 — |. 

Fig.4.a: tan/3 = 0.8, for three values of A 5 =4.0, 6.0 and 10.0. 
Fig.4.b: tan/3 = 1.6, for three values of A 5 =4.0, 8.0 and 12. 



Fig. 5 Fig.5.a: Bosonic contribution (T b ^ s ) to T# as function of A5. We chose: uih = 180, 
mh = 120, m H ± = 200, tjia = 110 (GeV), a = /3 — |, and several values of tan/3. 

Fig.5.b: Bosonic contribution (r^ s ) to as function of sin a for three different 
configurations Ci i=l, 2, 3: 

Ci : Wffi = 220, = 180, m/, = 80, tan/3 = 3.6 and A5 = 5 
C 2 : m H ± = 250, m H = 280, m h = 140, tan/3 = 1.6 and A 5 = 5 
C 3 : mH± = 420, m H = 400, m ft = 290, tan /3 = 2.6 and A 5 = 8 



Fig. 6 Ratio R (eq |5.1| ) as a function of tan f5 for various values of itla and for m#± = 500 
GeV. 
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